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ABSTRACT 


Measurements in the boundary layer and wake of a stalled airfoil are 
presented in two coordinate systems - one aligned with the airfoil chord, 
the other being conventional boundary layer coordinates. The NACA 4412 air- 
foil is studied at a single angle of attack corresponding to maximum lift, 
the Reynolds number based on chord being 1.5x10 . Turbulent boundary layer 
separation occurred at the 85% chord position. The two-dimensionality of 
the flow was documented and the momentum integral equation studied to illustrate 
the importance of turbulence contributions as separation is approached. The 
assumptions of simple eddy-viscosity and mixing- length turbulence models are 
checked directly against experiment. Curvature effects are found to be im- 
portant as separation is approached. 
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LIST OF SYMBOLS 


c 

c f 

C 

P 

H 

IX, IY 

P 

R 

R » 

u,v,w 

U 

e 

u * 
ref 

x,y 

k 

i 

t 

'Uj UJ 

oc 


s 

e 

v 

*^T 

f 


airfoil chord 

skin friction coefficient, *"C W / i p 
pressure coefficient, V> - h-* f / if 

boundary layer shape factor , 
indices referring to processed data mesh 
static pressure 

radius of curvature of airfoil surface 

Reynolds number based on momentum thickness, U J 

velocity components in (x,y) coordinates 

parallel component of velocity at edge of boundary layer 

reference velocity as measured by roof-mounted pitot-static tube 

coordinates normal and parallel to airfoil chord 

±( U-' M + V ,% -b 

mixing length, 'tf'v * 

K •+■ V'* -h U>' x 

velocity components in conventional boundary layer coordinates 0*,y. 
coordinates normal and parallel to local airfoil surface 
mean streamline direction relative to local airfoil surface, 

C v / ) 

boundary layer thickness: y * & when -u « O.^ 
displacement thickness: f f / — \ c/u 

f ~ J JT ) / 

momentum thickness: / ^ / / - u j ofv 

Ue ' UtJ * 

kinematic viscosity, f* 1 P 
eddy viscosity 
fluid density 
shear stress 
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LIST OF SYMBOLS (continued) 


Suffixes 

e 

w 


external stream conditions 
wall conditions 



1 . 0 INTRODUCTION 


One of the most detailed investigations of turbulent boundary layer separa- 
tion was that reported by Schubauer and Klebanoff (Reference 1). Conventional 
hot-wire anemometry limited their measurements to the region upstream of 
separation. Later, Bidwell (Reference 2) used this data to investigate the 
importance of turbulence terms in the momentum integral equation. 

Wadcock (Reference 3) made measurements in the boundary layer and near 
wake of a two-dimensional airfoil at a single angle of attack corresponding 
to maximum lift. The experiment employed a flying hot-wire to avoid direc- 
tional ambiguity in the hot-wire signals and so enable measurements to be 
made in the turbulent boundary layer downstream of separation. It is believed 
that these measurements will provide a major test case for development of 
calculation methods for turbulent flow, and so the above data (in the form 
of punched cards) have been filed with NASA-Ames Research Center. Details 
of the experimental configuration can be found in Reference 3. 

The data takes two essentially distinct forms: 

a) Raw Data 

The original output from the hot-wire inversion, consisting of data 
at discrete points along the flying-hot-wire arcs in a coordinate system based 
on the wind tunnel axis, was recorded on 14,620 punched cards (85 files, 

86 frames, 2 probes). A tabular description of the card format can be found 
in Reference 3. The raw data are thus defined at points which are awkardly 
placed. Further processing was therefore performed. 

b) Processed Data 

One objective was to use interpolation to redefine the data on a 
rectangular grid which is sufficiently fine to satisfy the analyst without 
driving the data beyond their real accuracy. Mostly for convenience in 
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describing the airfoil surface the grid in question is aligned with the air- 
foil chord. A further objective was to lay a foundation for the calculation 
of spatial derivatives. These final processed data were recorded on 29,195 
punched cards. Because the passage from raw to processed data was not only 
tedious but expensive, we expect most users will be willing to begin with the 
processed version of the data (see References 4 and 5) . Reference 5 contains 
a partial listing of the processed data. 

The rectangular grid used for the processed data is aligned with the air- 
foil chord. The method of interpolation and smoothing is detailed in Reference 
3. The mesh size is 1 cm in the chordwise direction and 0.2 cm in the cross- 
flow direction. Integer grid indices for the mesh, IX and IY , are related 
to the coordinates x and y of Figure 1 by 
IX = 1 + (x - 6.86) 

IY = 1 + 5 (y + 12.41) 

where x and y are in centimeters. IY ranges from 1 to 296 and IX from 1 to 
175 thereby documenting the wake flow for one chord length downstream of the 
airfoil trailing edge. 

2.0 PROCESSED DATA 

The processed data is presented in the form of contour plots in Figures 
2, 3 and 4. Figure 2 is a plot of the intermittency factor (fraction of the 
time that the flow is turbulent at a given position in space) in the boundary 
layer and near wake. The intermittency is close to unity over a substantial 
part of the boundary layer. To a lesser extent the same is true for the near 
wake. Figure 3 presents the mean velocity components normal and parallel to 
the airfoil chord. The reverse-flow region is clearly visible starting at 
x/c — 0.85 and extending to x/c — 1.07. The mean reverse-flow velocity in 
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the recirculation region is less than 207 o of the free-stream velocity. Figure 


4 presents data on the three double correlations u'u', u'v', and v'v*. In the 
boundary layer, as expected, velocity fluctuations parallel to the chord are 
very much larger than those normal to the wall. Downstream of the airfoil 
trailing edge, however, fluctuations in the two velocity components soon become 
comparable. The contour u'v' = 0 appears to emanate from close to the airfoil 
trailing edge, far from the apparent separation point as indicated by the 
mean velocity data. The explanation lies in the fact that at separation the 
shear stress is zero only when evaluated in boundary layer coordinates (i.e., 
normal and parallel to the local airfoil surface) . 

3 . 0 TWO-DIMENSIONALITY 

Before too much effort was expended in considering the success or failure 
of turbulence models in predicting the separated flow around the airfoil in 
question, it was necessary to document the two-dimensionality of the flow 
as fully as possible. Static pressure measurements on the airfoil surface 
(see Figure 5) and total pressure measurements through the wake approximately 
1% chords downstream of the airfoil trailing edge (see Figure 6) at three 
spanwise stations indicate excellent two-dimensionality. However, it is 
well known that pressure measurements are insensitive indicators of three- 
dimensionality- -a small lateral pressure gradient can result in substantial 
spanwise motion. 

The density of the data made available under the name "processed data" 
allows the computation of spatial gradients quite readily. It is instructive 
to test the two-dimensionality of the flow by examining the equation for 
conservation of mass: 

otuL _ 4 _ 

dK dZ 
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Clearly, if the left-hand side of the above equation is zero everywhere, 
then the flow is two-dimensional. If it is non-zero, then there will be cross- 
flow. 


Before computing s VaK and 3v /dy we need some measure of the reliability 
of the spatial gradients. This can be obtained by considering the vorticity 


%y _ du. 

a* ay 


3.1 VORTICITY 

In the external flow (i.e., outside both the boundary layer and wake) the 
flow is irrotational. The vorticity must be zero even if the flow is three- 
dimensional. This enables one to assign some sort of confidence level to 
spatial gradients of the mean velocity. Outside the boundary layer and wake 
such gradients are naturally much smaller than inside. However, useful infor- 
mation may be gathered from such an investigation. 

Working in the coordinate system based on the airfoil chord data is pre- 
sented in Figures 7, 8 and 9 along lines of constant x. 

(a) Through the boundary layer at IX = 50 (upstream of separation, 
x/c = 0.620) 

(b) Through the boundary layer at IX = 75 (downstream of separation, 
x/c = 0.897) 

(c) Through the wake at IX = 100 (x/c = 1.175). 

Spatial gradients are computed by simple differencing of the processed 
data (see Reference 3) which is available every 1 cm in x and every 0.2 cm 
in y. Just as all velocity components are normalized with the tunnel reference 
velocity U ^ , all spatial gradients are made non-dimensional by means of the 
airfoil chord, c. Because of the 5:1 ratio in grid step size in the arrays 
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containing the mean velocity data, gradients in the y-direction show five 
times the scatter displayed by the x-derivatives , i.e., x-derivatives appear 
much smoother. The absolute accuracy of the data is within 2 % for the dimen- 
sionless mean velocity components, but it is much harder to assign a similar 
figure to the spatial gradients. Examination of Figure 9 for x/c = 1.175 
shows that aa /»y bounces around somewhat on one side of the wake. This 
scatter is directly traceable to the raw data (see Reference 3) and unfortunately 
persisted throughout the smoothing and interpolation process. This problem 
is only present in the wake and does not appear in the boundary layer. The 
maximum scatter in the y-derivatives is therefore ^1. The x-derivatives are 
much smoother. From Figures 7, 8, and 9 it appears that the vorticity as 
calculated from is sufficiently close to zero in the external 

flow that the spatial gradients are reliable. 

3.2 CONSERVATION OF MASS 

Conservation of mass requires, for a general three-dimensional flow 

Su. ( o>V | 3w B o 

dK 3y dl 

We already have some degree of confidence in the spatial derivatives outside 
the boundary layer and wake. If we can assume that these gradients are equally 
reliable inside the boundary layer and wake (where they are necessarily much 
larger) then we can apply the above equation in the form 

+ «3>y — _ 2>w 

dK 9z 

to obtain some measure of the two-dimensionality of the flow (see Figures 10, 

11, and 12). Despite the difficulties inherent in the differentiation of 
experimental data, downstream of the airfoil trailing edge is zero within 
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reasonable limits considering the separately large values of and ®V 

At station 75 however, exhibits a departure typical of most boundary 

layer stations indicating some cross-flow or divergence. Absence of ex- 
perimental skin-friction data makes a momentum balance difficult to perform. 
Lack of an extensive logarithmic region in the velocity profile makes a 
Clauser plot for C^. unprofitable. It is therefore difficult to assess the 
degree of three-dimensionality with any certainty. 

4.0 DATA IN BOUNDARY LAYER COORDINATES 

Data in conventional boundary layer coordinates (normal and parallel to 
the local airfoil surface) were obtained approximately every 1 cm along the 
airfoil surface by interpolation in the processed data and rotation into the 
new coordinate system (see Figure 1). 

Figure 13 shows that the intermittency is close to unity over much of the 
boundary layer. Intermittency is a scalar quantity and hence is independent 
of the coordinate system chosen for presentation. Figure 13 and Figure 2 are 
therefore essentially the same. The mean velocity components tc and 1/ are 
shown in Figures 14 and 15. Figure 14 indicates that separation occurs close 
to station 70. The reverse-flow velocity is everywhere less than 20% of the 
free-stream velocity. Figure 16 presents the ratio of the mean velocity com- 
ponents ^ % c . It should be recalled that one of the Prandtl boundary layer 
assumptions is Perhaps a reasonable limit on the validity of such an 

assumption would be “** / vi = 0.1. As can be seen from Figure 16 such a value 
is reached well ahead of separation. Figure 17 presents the same data in a 
slightly different form. This figure shows the angle that the mean flow makes 
with the local airfoil surface. Note that for V / U. = 0.1, tan” 1 (^/it ) 

6 degrees so the flow does not have to depart from being parallel to the wall 
by very much to invalidate the assumption *0 ■« TC . 
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Figures 18, 19 and 20 present the three Reynolds stress terms tdli , 


tt'v' and 1/V*. As can be readily seen, downstream of separation very high 


turbulence levels are encountered with 




as high as 30%. Schubauer 


and Klebanoff (Reference 1) report much lower turbulence levels at separation, 
possibly due to poor instrument response (the time constants for their hot 
wires were approximately 0.002 seconds). Figure 19 shows that the corres- 
ponding velocity fluctuations in the normal direction are much smaller as 
expected. Figure 20 presents the Reynolds shear stress. The maximum shear 
stress is seen to occur first near the surface and move progressively out- 


wards with increasing downstream distance. At station 70 *U r , 'V > is close to 
zero near the wall as expected ( cf , Figures 3 and 4). However, starting slightly 


ahead of separation, V* becomes positive in the outer part of the shear 
layer. 

At first sight this was surprising, but a similar phenomenon has previously 
been observed with curved shear layers. Streamline curvature in the plane of 
the mean shear is known to produce surprisingly large changes in the turbulence 
structure of shear layers. To understand this effect, a short discussion is 
appropriate. A flow is said to be stable if, on displacement, a fluid particle 
is prevented from further displacement by a net restoring force. For flows 
over convex surfaces the centrifugal force is balanced by an inward pressure 
gradient. If a fluid particle is displaced outward across mean streamlines 
into regions of higher mean velocity, the centrifugal force on the fluid particle 
will be less than the new mean normal pressure gradient resulting in a re- 
storing force. Similarly, if the fluid particle is displaced inward, a net 
outward restoring force results. Thus convex boundary layer flows exert a 
stabilizing influence on turbulent momentum exchange resulting in a decrease 
in the interchange of both momentum and energy. 
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The direct effect of curvature on the Reynolds stresses was measured by 
So and Mellor (Reference 6) for a turbulent boundary layer on a convex wall. 

In their measurements the boundary layer was first developed on a flat surface 
and then passed to a convex surface; the R/£ ratio was about 12. The 
stabilizing effect of curvature was immediately evident resulting in the 
virtual disappearance of the turbulence in the outer layer. The ability of 
turbulence to produce shear stress was reduced by the stabilizing effect of 
curvature and a shear-stress profile that vanished in a region where the 
velocity gradient was still substantial was measured. A corresponding rapid 
decrease in skin-friction coefficient was also observed. 

A careful scrutiny of the data of So and Mellor at the most downstream 
stations reveals not simply the suppression of Reynolds shear stress but in- 
deed a sign reversal for in the outer half of the boundary layer. 

Their measurements proceeded some 24 boundary layer thicknesses downstream 
from the start of wall curvature in a flow of essentially zero pressure gradient 
with R/£ constant at about 12. It is not difficult to assume that if the 
measurements had continued downstream a more noticeable effect may have been 
visible. The present data is the result of the combined effects of stream- 
wise pressure gradient and continuously variable curvature parameters R/& . 

R /S is found to vary from 100 at the earliest stations f.or which boundary 
layer data is available to 20 at the airfoil trailing edge. The adverse 
external pressure gradient obscures the direct interpretation of curvature 
effects however. For example, the increased turbulent activity normally 
associated with the approach to separation apparently more than compensates 
for the stabilizing effect of curvature in the outer part of the boundary layer. 
No reduction in turbulent activity in the outer part of the layer is observed. 
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Thus we have a region in which s */3y and the Reynolds shear stress 
— are of opposite sign. The turbulent energy production term — t t'v' 
will have a negative value in this region. This implies that locally the 
energy of the mean flow is being increased at the expense of the turbulence. 

It is not the absolute values of the turbulent intensities but rather 
their ratio that provides useful information about changes in turbulence 
structure. Figure 21 presents the shear stress normalized with the turbulent 
energy and Figure 22 presents the same data normalized with the rms velocity 
fluctuations. ’U , V > f is a simple measure of the efficiency of the tur- 

bulence in producing shear stress. For a flat plate boundary layer its value 
is about 0.18 over most of the boundary layer decreasing rapidly near the edge 
of the layer. Measurements by So and Mellor (Reference 6) show that convex 
curvature lowers this value and, instead of going to zero at the edge of the 
velocity layer, it goes to zero at about half the boundary layer thickness. 

If we assume i £ 2 ^ ftV + then — * V A * * -ft) * v A*'* VvV i 

If the data presented in Figure 21 are examined in detail, we find that — it'v' J <£ 
is essentially constant across the boundary layer for the early stations. 

For example, at station 55 —’u'v /‘C takes the value 2/3 of 0.11 and at station 
60 the value 2/3 of 0.08, constant over the inner 75% of the layer. Down- 
stream of separation, of course, the character of this parameter is completely 
different due to the sign change associated with the Reynolds shear stress 
that occurs in the middle of the layer. 

For the flat plate boundary layer the correlation coefficient — 'M.'v ‘ V ' * 
is found to assume a constant value of 0.5 over the inner 3/4 of the boundary 
layer thickness, decreasing rapidly to zero at the edge of the layer. Convex 
curvature is found to reduce both the lateral extent of the region of constant 
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correlation coefficient and also the value of the coefficient (Reference 6) . 
Figure 22 shows that at the earliest boundary layer stations the correlation 
coefficient is already down to 0.4 and as separation is approached this value 
decreases rapidly. The assumption of constant correlation coefficient through 
the boundary layer also steadily becomes worse. 

Figure 23 presents data on the growth of the boundary layer thickness <S 
and includes information on the variation of the external velocity U. at 
the edge of the boundary layer. The boundary layer thickness increases rapidly 
as separation is approached, the growth rate assuming a constant value ^ Si O.J2 

dx. 

Figure 24 shows the development of the velocity profile every fifth station 
along the airfoil surface. At the earliest station for which boundary layer 
data is available (station 50, X’/c. * 0.620 ) a substantial velocity defect 
exists over much of the layer indicating the proximity of separation. Figure 
25 presents the velocity distribution through the boundary layer at each of 
the above stations the normalizing velocity being the local external velocity 
in each case. 

5 * At 

_ ^ __ , © and H = o /© , where 

the plane wall definitions were used in evaluating the integral thickness. 

c* 

Downstream of separation the displacement thickness e is seen to grow linearly 
with distance along the airfoil surface, x. — 0.24. At separation 

the shape factor H is approximately 3.6, Schubauer and Klebanoff quote a 

Jf 

shape factor close to 2.7 at separation. Table 1 lists o , © and H at a 
few select stations. No data is at present available on the surface shear 
stress distribution. 
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5.0 THE REYNOLDS EQUATIONS 


The motivation behind presentation of data in boundary layer coordinates 
was the term-by-term examination of the time-mean Navier Stokes equations in 
order to investigate the validity of the Prandtl boundary-layer assumptions 
in the vicinity of separation. 

The Navier Stokes equations for a two-dimensional flow in which the vis- 
cosity and density are assumed constant and the body forces are negligible 
are : 


U«5ul V du 

at ax 5y 

■a - 1 

7 

e>^ 4 . ^ ( a a u w ^ 

5 x ' 5 k v ay*/ 

(i) 

3v -v v dv 

at ax 5y 

u 

i 

■tl- 

aw 4 . v ( a x v a*y \ 

ay ' ax* ay 1 / 

(2) 

The associated continuity equation is 



3 u. 

ax 

-v 

a y 

— 0 

(3) 


By resolving the instantaneous values into mean and fluctuating components, 
by considering the mean flow steady and taking a long-time average of the 
fluctuating components and by using the continuity equation, equations (1) 
and (2) may be written as 


u. -v V a u ) a _q. V / i'ul d*uL 

ax *y ax ay f Ax V AX* dy*/ 


c: ay v ay + a ^u'v 1 \ a vV 1 - - i at + ■»> / a l v a 2 v \ (5) 

ax ay ax ay f Sy \ ax* ay 1 / 
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Considering the Prandtl boundary- layer assumptions valid (d « d v <« iuL ), 

an iy 


vL du 

*x 


« + a lu‘u') + a ^‘v’) * -! at 
ay ax ay f 5> x »y* 


( 6 ) 




f 4 


( 7 ) 


Combining the above two equations. 


VX du V ac. a TxV — vV J c Ug, cUde. * 4 - 1 

ax ay ax «=a* f ay 


(8) 


where U e is the velocity at the edge of the boundary layer and '“C * «V 'TT-j 

^ = - f ^ . 

The above equations were written for flow along a flat wall. For two- 
dimensional flow along a curved wall we can write the corresponding equations 
in a curvilinear orthogonal system of co-ordinates whose .^c-axis is parallel 
to the wall, the t^-axis being perpendicular to it (see Figure 1). With the 
assumptions that the boundary layer thickness is small compared with the 
radius of curvature of the wall and viscous terms are small compared with 
inertial terms (i,e., high Reynolds Number) we obtain 


Bu + VL 9tt v du. ILI/ c 

- 1 

(la) 

&k. Bx. Sy 

f 3c 



By + uBy c - j L BJ> 

<?£. Bx. ft f 


(2a) 
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0 >U 

S>x. 


(3a) 


+ y. » o 

For the flow under consideration the curvature terms are important. The 
object of the following discussion is to establish the importance of individual 
terms in the time-mean version of the above equations , namely 


as>u. + 'V ss + & ( •u'uy + & ( * t v') + /£SL±z±!*l) 

dx* &X. <fy ' Jp ' 


-/ dt 


(4a) 


u 

<2* 


v&y +. £ ( zt'Vj +. f'V f 'V , \ - f Hu + 'u , u > - yV ; ) _ dp 


(5a) 


Each term on the left-hand side of equations (4a) and (5a) has been 
evaluated at three streamwise positions on the airfoil surface: 

Station Description 

55 Upstream of separation 

70 Close to separation 

80 Downstream of separation 

Spatial gradients were estimated by simple differencing in data stored at 
grid locations with a mesh size of approximately 1 cm by 0.2 cm in the ortho- 
gonal directions and y , parallel and normal to the local airfoil surface 
respectively. 

Figures 27, 28 and 29 present a study of individual terms in the y— 
momentum equation (equation 5a) at each of the above three stations x/c = 0.675, 


17 



0.842 and 0.952. It is clear from all three figures that curvature terms are 


important at all stations along the surface of the wing. The usual neglect 


of the terms and BBv from the y-raomentum equation is seen to have little 

justification. Each term is separately quite large and comparable with the 

term that is retained, namely . The one term that can quite safely 

y 

be neglected throughout the boundary layer is $ (h’v'} . 

Sx 

At the outside edge of the boundary layer ?r t£ O for stations 55 and 70 

7 


whereas downstream of separation (station 80) this is not the case. Neglecting 
turbulent contributions to the momentum equation, equation (5a) becomes 


or 


“U By 



•v Bv — -u a - — J Bp 

7Z P 



Large transverse pressure gradients are thus to be expected in regions of rapid 
streamline curvature, e.g. close to separation where the mean streamlines are 
closed. However, in the recirculation region the mean velocity is low and 
hence large transverse pressure gradients are not encountered. The quantity 
in the square brackets in the above equation represents the difference be- 
tween the mean streamline curvature and the local curvature of the wall. 

Thus we see from Figures 27 and 28 that the external flow is moving parallel 
to the wall at stations 55 and 70. Downstream of separation however, Figure 
29 indicates that the external flow is no longer parallel to the local airfoil 
surface indicating that the shear layer has simply peeled off. Examination 
of Figures 27, 28 and 29 shows that the largest pressure difference across 
the layer for the three stations considered occurs at station 80. This is 
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the result of having both the largest transverse pressure gradient and the 


thickest layer present at this station. Integration of across the layer 
at this station yields a value for ^ = <V C L of approximately 0.07. 

Figures 30 , 31 and 32 refer to individual terms in the ^.-momentum equa- 
tion (equation 4a) evaluated at stations 55, 70 and 80 corresponding to x/c = 
0.675, 0.842 and 0.952. Figure 30 shows that ahead of separation (station 55) 
a large adverse pressure gradient is present throughout the boundary layer. 
Figure 32 shows that at station 80, downstream of separation, °P > © in the outer 

Sk 

half of the layer whereas close to the wall dp < © . This is entirely con- 


sistent with the fact that fluid close to the wall is moving upstream at this 
station. As fluid moves upstream it must move into a region of increasing 
pressure in order to be brought to rest at the separation point. Thus there 
should be a favorable streamwise pressure gradient close to the wall at station 
80. The streamwise pressure gradient at the wall has also been evaluated from 
the surface pressure distribution shown in Figure 5. This information is also 
shown plotted in Figures 30, 31 and 32 for comparison. Note that at the wall, 
equation (4a) reduces to 


- L ) = 2 / 'u , 'v' ) 

r **/„ ^ 

so that comparison between the pressure gradient evaluated from the flying- 
hot-wire measurements using equation (4a) and that computed from surface 
pressure measurements is really only a comparison between the streamwise 
gradient of the surface pressure and the normal gradient of the Reynolds 
shear stress. 
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6.0 THE MOMENTUM INTEGRAL EQUATION 


Figure 33 examines individual terms in the momentum integral equation, 
neglecting contributions from turbulent fluctuations: 

3 * = 4? + ( H + 2) e dUe. 

c/jc. Ue. dx. 

If the above equation is used to describe the variation of the wall shear 
stress we move downstream and approach separation rc w is found to in- 

crease (at least initially) . This is clearly incorrect as the boundary layer 
is growing in an adverse pressure gradient (see Figure 23 for U e C'-)) and so 
should fall monotonically as separation is approached. This problem has been 
encountered by other researchers (see reference 2). Downstream of station 

72 ' T iv < ° indicating reversed flow cf mean velocity data shown in Figure 

14. The agreement appears fortuitous considering the indicated behavior of 
‘■'W ahead of separation. 

Following Bidwell (Reference 2) , the momentum integral equation can be 
written 


'*T W = dS + (H+2) 9 dM g 
pul dx Ue dx. 



L 





2 C *V; + £> (v'*') ) c/ydu 

** / t 4 



where curvature effects have been neglected and only turbulent contributions 
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to “&P retained in equation (5a) . Contributions to the momentum integral 
equation from the above turbulence terms were evaluated at a limited number 
of stations. The results are presented in Table 2. Clearly the most signifi- 
cant contribution comes from the streamwise derivative of . This term 
only involves the first derivative of and hence can be determined relatively 

accurately from the experimental data. Contributions from other terms are 
much smaller, although still significant, and rely on the calculation of second 

derivatives their accuracy is therefore diminished. The cumulative effect 

of these turbulent contributions is insufficient to account for the growth of 
Assuming that the mean flow is indeed two-dimensional this can only 
be due to the omission of curvature terms from the equations of motion and the 
neglect of mean convective terms in the ^-momentum equation. 

7.0 TURBULENCE MODELS 

There are two essential problems associated with Reynolds averaging 
(Reference 7). First, since the Reynolds stresses are unknowns, Reynolds 
averaging leads to a set of equations with more variables than equations. An 
attempt to obtain equations for the. Reynolds stresses by successive multipli- 
cation of the Navier Stokes equations followed by time averaging leads to a 
hierarchy of equations in which, for example, 'U / V / depends on terms of the 
form 'u’xt'v* and so on. The hope in using Reynolds averaging lies with the 
possibility of closing the resulting sequence of equations by a physical or 
at least a statistical argument. This is the so-called closure problem. 

The second problem with Reynolds averaging is related to the usefulness 
of such averaging. Recently it has become clear that Reynolds averaging is 
not as obviously useful as has previously been thought. Time-resolved measurements 
of the statistical properties of such terms as l**V # show, at least in some 
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cases, a nearly bi-modal probability distribution; the dimensionless coefficient 
— iaJv* jJ is found to be zero for stretches of time interrupted by 
bursts in which it is unity; its mean value, the correlation coefficient, is 
of the order of 0.5. The physical significance of the mean value of a variable 
fluctuating in this fashion is at least questionable. 

If the viscous shear and normal turbulent stress terms are neglected in 
equation (8) and the resulting equation combined with the time-mean of equa- 
tion (3) we have a system of equations that can be solved for any given ex- 
ternal pressure gradient, provided something can be said about the turbulent 
shear stress — ^ t/v' . The primary objective of all turbulence models is 
therefore a description of the Reynolds shear stress. 

The simplest models relate the turbulent shear stress uniquely to the 
mean flow conditions at each point. Since these models require only algebraic 
expressions they have enjoyed prolonged popularity. All models of this class 
use the eddy -viscosity concept of Boussinesq: 

— n { v' sjy- 3a 


where V-y- is the eddy viscosity. The two most familiar examples of this 
class are distinguished by the way in which is calculated. 

a) Cons tant-Eddy -Viscosity Model, =r C. ^ C ^«ax " 

All quantities on the right-hand side of this equation, apart from , 
may be functions of the longitudinal distance X ; '3^- is supposed uniform over 
any cross-section. 


b) Mixing Length Hypothesi 


Ls, . 


/ * x 

JL is usually taken as proportional to o . 


The mixing length JL is usually taken as proportional to d . This 
method is still the basis of many calculations of the turbulent boundary layer 
which are carried out today. 
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For many boundary layer flows, Prandtl's mixing length hypothesis works 
surprisingly well. The spreading rate, as well as velocity profile, can be 
predicted satisfactorily. However, the constants involved must vary with the 
problem considered. This lack of universality is an indication that the under- 
lying model of turbulence lacks some important features of real flows. 

Figure 34 plots the ratio — against Vs for four stations ahead 

of separation. Clearly, the closer the approach to separation, the less 
accurate the assumption - Con&hanb becomes. The constant-eddy-viscosity 
model assumes that scales with S and Ue. ; for example, where 

S * S (x.) , Ue C*) and Ot is a constant. Figure 35 shows that there is 

a substantial decrease in as separation is approached. Galbraith and 

Head (Reference 8) found the same to be true for the experimental data of 
Schubauer and Spangenberg , Flow B (Reference 9) . They found (iS T / ) MAKt 
decreased from approximately 0.018 to 0.007 as separation was approached. The 
present data apparently starts where the Schubauer and Spangenberg data finishes 
since ) is already down to 0.007 at the earliest station for which 

' 'max. 

boundary layer data is available. Various flows were examined by Galbraith 

and Head revealing values for as high as 0.028. A particular 

/ , * ' MAK 

value for L+rU**) may be typical for a given flow but is apparently 

MAfL 

highly dependent on the state of equilibrium of the flow. 

As mentioned above. Figure 34 shows that as separation is approached, 
the assumption of constant eddy viscosity across the boundary layer becomes 
less accurate. Examination of Figures 14 and 20 reveals that downstream of 
separation the constant-eddy-viscosity model fails completely. Whereas 
is positive throughout the boundary layer (apart from a very thin layer close 
to the wall) — 14* V * is positive close to the wall and negative in the outer 
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half of the layer. The ratio -'u'v' f <du changes sign across the boundary layer 

/ 

and is thus far from being constant. It therefore appears that the constant- 
eddy-viscosity model has limited success well ahead of separation, but as 
separation is approached the model is unrealistic. 

If we examine the mixing length hypothesis the story is somewhat similar. 

The mixing length is given by Jl - -t / . Figure 36 shows that at 

• < 9 ^ / * 5 ^/ 

the earliest stations the assumption 1 = cons ban b may not be too bad, but 
the assumption gets progressively worse as separation is approached. The 
mixing length is assumed to scale with <£ for the outer part of the boundary 
layer. Figure 37 shows that this assumption is valid well ahead of separation. 
Close to the wall the mixing length is generally assumed to be described by 
A The agreement is satisfactory although insufficient data 

is available close to the wall to be conclusive. For many boundary layer flows 
S - O'/C for S £ G'2.15 (see Reference 8). However, for the early boundary 
layer stations a typical value for is 0.05, half the above value. Gal- 

braith and Head (Reference 8) also noticed a decrease in ■i/s for the outer 
part of the layer as separation is approached. Downstream of separation the 
mixing length model fails for essentially the same reason the constant-eddy- 
viscosity model failed. Once again, in boundary layer coordinates, is 

positive throughout the layer (except very close to the wall) whereas 
changes sign in the middle of the layer. A negative value for the ratio 
is not consistent with the physical model. 


’u'v' / jDu I i 

/ I I 

One of the more successful turbulence models employs the turbulent kinetic 


energy equation and assumes that the shear stress < Zj‘- - is closely re- 
lated to the turbulence kinetic energy £ £ I*'* V -h LO ,X ^ Townsend 

proposed that the ratio xl'v* to M - i (*'*->■ ** + <*") might be a universal 
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constant. Bradshaw employed this assumption with great success as part of a 
calculation procedure for attached boundary layers. For many flows, however, 
the underlying assumption is not realistic. Often the shear stress vanishes 
where the turbulence kinetic energy iz remains finite, or, as in the case 
for round jets, where k. has its maximum value. Bradshaw defines the Reynolds 
shear stress by --u'v's 4 <^ 2 where c£ * ( V* *+* V'V co'y and assigns the 


value 0.15 to CL. 


Examination of Figure 21 shows that this is reasonable 


only for the very earliest of boundary layer stations. As separation is 
approached the assumption becomes less accurate. 

Yet another model includes the proposal — zi't/' * > 7 - where the eddy 

% L ^ 

viscosity is allowed to vary through the layer according to s C M L. 
and jk = l •+ 40 'M is calculated from a differential equation and 

the length scale from an empirical algebraic expression. This model suffers 
from the same problem the constant-eddy -viscosity model has namely, down- 

stream of separation kV' passes through zero in the middle of the boundary 
layer where is far from zero. 


It should be emphasized at this point that the major problem with all the 
above turbulence models is associated with the change in sign of the Reynolds 
shear stress in the outer part of the boundary layer downstream of separation. 
If the hot-wire data is examined in the original coordinate system of Figure 

4 (axes normal and parallel to the airfoil chord) , this problem no longer 

/ He 

u t s is more nearly constant across the layer, 
even downstream of separation, although the problem of predicting the value 
of 0 C remains unsolved. Perhaps this is an indication that a turbulence 
model may fail because of an inappropriate choice of coordinate system. 

Close to separation, boundary layer coordinates are perhaps inappropriate 
anyway. An attractive alternative coordinate system that merits attention is 
orthogonal to the local mean streamlines. 
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Bradshaw (Reference 10) discusses the effect of "extra" strain rates 


(i.e., other than a simple shear ). For a small extra strain rate e, 

*9 „ « 

the Reynolds stresses are expected to change by a factor F = 1 + & . — 

> ' 

where f? varies from case to case but is of order 10. The rate of strain 

is defined as small if e-/ (_3u/&j) is numerically less than 0.05 so that 
0.5 < F < 1.5. Outside this range a linear correction factor like F cannot 


be trusted. For a thin curved shear layer the extra rate of strain is 6. 
in (x,y) coordinates or e = -U/(y+A) in coordinates. To account 

for the observed change in sign of the Reynolds shear stress requires F to 
change sign, which takes us outside the range of validity of the model. 

8.0 CONCLUSIONS 


Sx 


A shortage of experimental data limits the development of turbulence 
models. The determination of the most economical computational model for the 
prediction of subsonic flows separating from solid surfaces remains a problem 
of considerable interest. Numerical solution of the full Navier Stokes equa- 
tions requires considerable computer time and thus provides motivation for 
considering the range of applicability of methods based on more approximate 
forms of the governing equations. 

The richness of the experimental data was illustrated by the use of 
contour plots. The reliability of spatial gradients was also documented. 
Individual terms in the time-mean Navier Stokes equations were examined to 
determine their relative importance. In particular, the usual simplification 
of the y-raomentura equation to s - ^v'* was shown to be unjustifiable close 
to separation. 

Contributions to the momentum integral equation from turbulence terms 
were shown to be essential as separation is approached. The form of the 
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momentum integral equation used by Bidwell failed to represent the behavior 
of the wall shear stress even with turbulence contributions included. This 
was attributed to the neglect of certain terms from the y-momentum equation, 
terms which are large close to separation. 

Eddy-viscosity and mixing-length turbulence models are closely-related 
hypotheses of a connection between the Reynolds shear stress and the local 
mean velocity gradient. Mixing length theory with its concept of momentum 
transfer by fairly well organized lumps of fluid is considered a good first 
approximation to the behavior of turbulence dominated by large eddies. The 
large eddies carry most of the shear stress and their lifetime is so long 
that during it they may travel downstream for a considerable distance. Thus 
the Reynolds stress at a given position depends significantly on upstream 
history and is not, for example, uniquely determined by the local mean velocity 
gradient as in a laminar flow. This illustrates the need for transport 
equations to describe turbulence. The main reason for the success of mixing- 
length and eddy-viscosity is that many fluid flows used as test cases for 
turbulence models are nearly in a state of self-preservation or local equilib- 
rium where profiles of mean velocity, Reynolds stresses, etc. are geometrically 
similar. Under such conditions both models work well. Complex flows, far 
from equilibrium, fare less well. In particular, for the airfoil boundary 
layer under consideration, mixing- length and eddy -viscosity models fail to 
predict levels of *'£’ and cannot handle the outer part of the boundary layer 
where curvature effects appear to dominate. Clearly, more sophisticated 
turbulence models are required in order to predict such a complex flow. 
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TABLE 1 


BOUNDARY LAYER INTEGRAL PARAMETERS 


Station 

x/ c 

,* 

100 < /c 

100 ©/c 

H 

V U ref 

*0 

50 

0.620 

0.827 

0.393 

2.10 

1.295 

7750 

55 

0.675 

1.093 

0.479 

2.28 

1.247 

9110 

60 

0.731 

1.526 

0.600 

2.54 

1.207 

11040 

65 

0.786 

2.230 

0.756 

2.95 

1.178 

13560 

70 

0.842 

3.288 

0.913 

3.60 

1.158 

16110 

75 

0.897 

4.536 

0.994 

4.56 

1.134 

17170 

80 

0.953 

5.845 

1.014 

5.76 

1.112 

17180 


c = 0.9012 m 

U . = 27.13 m/s 
ref 

^ = 0.1605x10 ^ m~ / s 

R © = D e®/ V 
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TABLE 2 


TURBULENCE CONTRIBUTIONS TO THE MOMENTUM INTEGRAL EQUATION 


Station 

x/ c 


A 

B 

C 

55 

0.675 


-0.004 

0.000 

0.001 

60 

0.731 


-0.006 

0.000 

0.001 

65 

0.786 


-0.013 

-0.001 

0.003 

70 

0.842 


-0.015 

0.001 

0.001 

75 

0.897 


-0.012 

0.001 

0.002 

80 

0.953 


-0.011 

-0.003 

0.006 



A: 
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NACA 4412 AIRFOIL SECTION 


Figure 1. Airfoil geometry and orientation of axes. 
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Figure 4. Contour plots of dimensionless double corr elations from final 
processed data. Contour intervals: u'u J ; o.oio 

u' v’ j ; o.oo4 

T7/u^ o.o io 





FLOW GUIDE 


/ 



(a) 


H/Qref 



Figure 6. (a) Configuration of model and wake. 

View is inverted from actual 
orientation in tunnel. 

(b) Two -dimensionality of wake with 
flow guides in optimum position. 

H is total pressure in wake referred 
to atmospheric pressure. 

Measurement position is 1.28 chords 
downstream of airfoil trailing edge. 
Airfoil is lifting in direction of 
increasing Y. 
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(b) 
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Figure 7 
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Figure 8. Vorticity distribution at IX = 75 , x/c = 0.897; 


3iy 
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Figure 10. Conservation of mass at IX = 50, x/c - 0.620; 


du 

ax 


< 3v 


-3w . 
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DISTANCE NORMAL TO AIRFOIL CHORD 


Figure 12. Conservation of mass at IX = 100, x/c = 1.175; 


DISTANCE NORMAL TO VW 



STATION 


Figure 13. Contour plot of intermittency factor. Contour interval is 0.1 



STATION 


Figure 14. 


Contour plot of mean-velocity component parallel to surface. 
Contour interval 0.1. 
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Figure 15. Contour plot of mean-velocity component normal to surface. 
Contour interval 0.1. 
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Figure 16. 


Contour plot of mean-velocity component ratio. 
Contour interval 0.1. 
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Figure 17. Contour plot of mean flow direction, 0C° 
Contour interval 2 degrees. 


bau," (V/U ). 
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STATION 


Figure 18. Contour plot of double correlation, iaJvl' f 
Contour interval 0.01. 
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STATION 


Figure 19. Contour plot of double correlation 
Contour interval 0.01. 




Figure 20. Contour plot of double correlation, 'LL V* f . 

Contour interval 0.001. 
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Figure 21. 


Contour plot of correlation coefficient, xt' l/‘ j '( Zl’ it' -*■ i/’v'J 
Contour interval 0.1. 
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Figure 22. 


Contour plot of correlation coefficient, 
Contour interval 0.1. 
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Figure 26. Variation of displacement thickness, momentum thickness 
and shape factor along airfoil surface. 
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Figure 29. Examination of time-mean y -momentum equation 
at Station 80. See Figure 27 for legend. 
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Figure 32. Examination of time-mean 3C- -momentum equation 
at Station 80. See Figure 30 for .legend. 
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Figure 33. Examination of momentum integral equation 
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Figure 37. Variation of mixing length through attached boundary layer. 
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